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We consider a non-spherical colloidal particle immersed in a fluid close to its critical point. The 
temperature dependence of the corresponding order parameter profile is calculated explicitly. We 
perform a systematic expansion of the order parameter profile in powers of the local curvatures of the 
surface of the colloidal particle. This curvature expansion reduces to the short distance expansion 
of the order parameter profile in the case that the solvent is at the critical composition. 
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I. INTRODUCTION 

Various interactions between colloidal particles im- 
mersed in a solvent can influence the properties of the 
system. In a colloidal suspension containing particles of 
different size, there is an effective attractive interaction 
between larger particles. This attraction is due to the 
extra volume that becomes available to smaller parti- 
cles when larger particles approach each other, leading 
to an overlap of excluded volumes which increases the 
entropy of the system (see, e.g., Ref. [l| and references 
therein). Colloidal particles, surfaces of which dissociate 
in solution, exhibit a screened Coulomb interaction due 
to charged surfaces of the particles and the surrounding 
counterions (see, e.g., Ref. Van der Waals disper- 

sion forces arise from induced dipole-dipole interactions 
due to quantum mechanical fluctuations of the charge 
density (see, e.g., Ref. [sj). The confinement of critical 
order parameter fluctuations in a binary liquid mixture 
near its critical demixing point give rise to long-ranged 
critical Casimir forces between immersed colloidal par- 
ticles (see, e.g., Ref. [1]). The richness of the physi- 
cal properties of colloidal suspensions is mainly based 
on the possibility to tune these interactions which dif- 
fer significantly in strength and range. In the case of 
the effective entropic interaction, the Coulomb interac- 
tion, and the van der Waals interaction this tuning is 
accomplished by changing the composition of the sol- 
vent by adding depletion agents, salt, or other compo- 
nents. For example, by matching the indices of refrac- 
tion of the colloidal particles and the solvent it is pos- 
sible to effectively switch off the dispersion forces and 
to create colloidal suspensions in which the actual ef- 
fective interaction between uncharged colloidal particles 
very closely resembles a hard core potential. In the case 
of charged colloidal particles, even small polyelectrolyte 
additives can have substantial impact on the aggregation 
and kinetic stability of the charged particles |5j]. Com- 
pared with such modifications, changes of the tempera- 
ture or pressure typically result only in minor changes 
of the effective interaction between colloidal particles. 
However, effective interactions generated by bringing sol- 
vents close to a phase transition of their own are ex- 



tremely sensitive to such changes as observed experimen- 
tally 1, 0, i, S M El, m, m, lii, m, m . in these exper- 
iments, a dilute suspension of spherical colloidal particles 
has been formed within the one-phase region of a binary 
liquid mixture acting as the solvent. As the temperature 
of the system approaches the phase separation tempera- 
ture of the binary liquid mixture, the colloidal particles 
aggregate and flocculate out of solution. Light scatter- 
ing measurements indicate that an adsorption layer rich 
in one of the two solvent species forms around each col- 
loidal particle. The experimental results strongly suggest 
that the colloidal aggregation behavior is induced by the 
presence of this adsorption layer. The thickness of the 
layer increases upon approaching the aggregation line in 
the phase diagram. 

In a classical binary liquid mixture near its critical 
demixing point, the order parameter is a suitable con- 
centration difference between the two species forming the 
liquid. The generic preference of confining boundaries 
for one of the two species results in the presence of effec- 
tive surface fields leading to nonvanishing order param- 
eter profiles even in the one-phase region of the phase 
diagram [13, These critical adsorption profiles be- 
come particularly long-ranged due to correlation effects 
induced by the critical fiuctuations of the order param- 
eter of the solvent. While several theoretical investiga- 
tions have been devoted to the understanding of critical 
adsorption phenomena on planar walls and spherical par- 
ticles (see Refs. [1, [l^ and references therein), critical 
adsorption on non-spherical particles has been studied 
only in the limiting cases of infinitely long cylinders (20| 
and very small dumbbells and ellipsoids within the frame- 
work of a small particle operator expansion [2l| . despite 
the growing interest in non-spherical colloidal particles 
(see, e.g., Refs. [HHl])- to the reduced symmetry of 
the shapes of non-spherical colloidal particles compared 
with spherical symmetry the interactions between them 
depend not only on the separation between their cen- 
ters but also on their mutual orientations. As a result 
various properties of fiuids consisting of such particles 
differ from the corresponding ones of fluids consisting of 
spherical particles. As a prerequisite for studies - mo- 
tivated by the aforementioned aggregation phenomena 



2 



near criticality ~ on critical Casimir forces between non- 
spherical colloidal particles and surfaces we investigate 
here systematically the temperature dependence of crit- 
ical adsorption on a single colloidal particle. In order to 
treat ellipsoids, spheres, and cylinders in a unified way 
within the appropriate field-theoretical approach and for 
general embedding dimensions D it is helpful to consider 
the particle shape of a hypercylinder defined as 



/Cd({i?J) - 




where Ri < R2 < ■ ■ ■ < Rd are the semi-axes of the 
particle (see Fig. [T]). In the case of equal semi-axes Ri — 
R2 = ■ ■ ■ = Rd the hypercylinder reduces to the so called 
generalized cylinder with an infinitely extended "axis" 
of dimension D — d [13, . In D = 3 the " axis" can 
be the axis of an ordinary infinitely elongated cylinder 
{d = 2), or the midplane of slab {d = 1), or the center 
of a sphere (d = 3). In the case that not all semi-axes 
are equal but d = D the hypercylinder reduces to an 
ellipsoid which is called a spheroid if the lengths of two 
semi-axes are the same. In Z) = 3 one may distinguish 
prolate spheroids (i?i = i?2 < R3) from oblate spheroids 
(i?i < i?2 = Rs). The general elhpsoid (-Ri < i?2 < 
R3) is called a triaxial ellipsoid. The generalization of D 
to values different from three is introduced for technical 
reasons because Due = 4 is the upper critical dimension 
for the relevance of fluctuations of the order parameter 
leading to a behavior different from that obtained from 
mean field theory valid in Z) = 4. It proves convenient 
to express the position vector r in terms of a distance 
r±_ perpendicular to the surface and dimensionless angles 
{Os}, s ~ 2, - ■ ■ ,D (see Fig.[T]). Moreover, the shape of a 
hypercylinder can be characterized by the smallest semi- 
axis R=Ri and the dimensionless ratios {6s = Rs/R}, 
s = 2, • • • ,d. 

II. ORDER PARAMETER PROFILES 

A. Scaling properties and short distance expansion 

Close to Tc the critical adsorption on the surface of 
the mesoscopic particle is characterized by an order pa- 
rameter profile {4>{r))t which takes the following scaling 
form 

Wr))t = a\tfP±irJ^±,R/(±,{9s},{5s}) (2) 

for distances r± from the surface sufficiently larger than a 
typical microscopic length. Here (• • • )t denotes the ther- 
mal average. ^±{t ^ 0) = Co^l^l is the bulk correlation 
length above (+) or below (— ) the critical temperature 
Tc, where t = [T —Tc) /T^ is the reduced temperature and 




FIG. 1: Schematic side view of an ellipsoid with semi-axes 
R = Ri < R2 < ■ ■ ■ < Rd- Only the projection of the particle 
onto the xi — Xd plane is shown. Any point r outside the 
particle can be reached by a distance r± measured from the 
point closest to r on the surface of the particle. The angle 
between the position vector of the latter surface point and 
the Xd axis is denoted as 9 = Bd- 



(3 and u are the standard bulk critical exponents. The 
scaling functions P±{r±/^±,R/£,±, {Og}, {Ss}) are univer- 
sal once the nonuniversal bulk amplitudes a and are 
fixed by the value {(j)b)t^o- — a\t\^ of the order pa- 
rameter in the unbounded bulk and by the true corre- 
lation length defined by the exponential decay of the 
bulk two-point correlation function in real space. There- 
fore on finds P_(r_L/C- ^ 00, R/^-, {9^}, {Ss}) = 1 and 
P+{r±/£.+ 00, {Os}, {Ss}) = 0. In the opposite 

limit r±/£_± 0, i.e., T — s- Tc, the scaling functions 
and the order parameter profile exhibit short-distance 
singularities in the form of power laws which refiect the 
anomalous scaling dimension of the order parameter: 

(^(r)),.o = «C±(r^/i?, {Os}, {Ss}) (^^ j . (3) 

The ratio P/v of the critical exponents has the value 1 in 
Ti = 4, « 0.517 in T) = 3 and 1/8 in £> = 2. The am- 
plitude functions C'±{r±/ R^ {^s jj {Ss}) are universal but 
depend on the definition of the correlation length, as the 
scaling functions P±{r±/^±,R/^±,{6s},{Ss}), too. Us- 
ing the operator-product expansion (see, e.g., Refs. f2^, 
|20|) a short distance expansion of the universal amplitude 
functions can be derived: 

C±ir^/R^O,{0s}ASs}) 

c± R 

+ [XkK + Xh2H'] (^)' 

+ [Xh^H^ + XhkHK + XaG] (^)' + O (^)' , 

(4) 

where H = iZ({0 J, {(5,}), K = K {{9 s} , {S s}) , 
and G = G{{9s},{5s}) are dimensionless local 
curvatures characterizing the surface and A/ with 
/ = H,K,H^,H^,HK,G are dimensionless coefficients 
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which depend on D but not on the shape of the bound- 
ary surface, i.e., they are independent of {Og} and {Sg}- 
The dimensionless local curvatures are related to the lo- 
cal radii of curvature Ai = Ai{{9s}, {5s}) according to 
(see, e.g., Ref. [28]) 
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(5b) 



i<j<k 



In the limit r±/R the universal amplitude func- 
tions C±{r±_/R = 0, {6*5}, {(5s}) reduce to the uni- 
versal amplitudes c± for the critical adsorption pro- 
files near planar walls confining semi-infinite systems 
[13. According to Eq. © C_ (r^/i?, {6IJ, {^J) = 
{^tli^Y'''C+{r^/RA0s}A5s}), where H.Y'^ = 
C- / c-\. « 1.4 (for D = 3) is a universal amplitude ra- 
tio which depends on the embedding dimension D (see 
Table I in Ref. [131 )■ Therefore it is sufficient to consider 
only one of these functions. 

In the case of a sphere, i.e., iJ^'P'*' ^ D ~ I, K'^'p''^ = 
{D-l){D-2)/2, and G'^'^p''^ = (L»- l)(£'-2)(D-3)/6, 
the universal amplitude functions are known exactly for 
any spatial dimension D by means of a finite conformal 
mapping from the half-space [29| : 



d^^"\r^/R) = c+ 



1 



2 R 



(6) 



A comparison of the short distance expansion of 
C^^P^^ir^^jR) with Eq. Q yields 

1/3 



2 V 



(7) 



and m. 



as well as 
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In D — 2 all coefficients A/ with I = H, , can 
be deduced from Eqs. ([7]) - (O. For D > 2 & determi- 
nation of the coefficients appearing in Eqs. ^ and Q 
would require the additional knowledge of universal am- 
plitude functions around differently shaped walls. How- 
ever, these functions are not available beyond the mean 
field approach discussed in Sec. IIIII below. 



B. Curvature expansion 

Recently a curvature expansion of the density profile of 
a hard sphere fiuid close to a big convex particle has been 
proposed and successfully applied to the density profile 
around a hard ellipsoid 30]. The curvature expansion 
separates the properties of the fluid from the geometry 
of the big particle and allows one to determine density 
profiles in complex geometries based on those obtained 
in much simpler geometric configurations. Assuming that 
an analogous curvature expansion holds for the order pa- 
rameter profile ((/)(r))t leads to 

{<t>{v))t^<t>'f\r^) + <j>f\r^,R)H 
+ 0f ) (r^ , R)K4>\"' ) (r^ ,R)H^ + <t>["' ^ [r^ , R)H^ 
+ {r±,R)HK + 4^\r^,R)G + 0{R-^), (10) 

where (j)f'^ (rx) is the order parameter profile for the half- 
space bounded by a planar wall and (r^ , R) with / = 
K, H^, HK, H^, G are expansion coefficient functions. 
A comparison of this equation with Eqs. ([3]) and ^ leads 
to the following short distance behavior of the curvature 
expansion coefficient functions: 



b\'\r^,R) 



-ac± 



Co" 



(11) 



-ac± 



■A,^i'^(rx/e±), 



(0) 



1, 



(12) 



for 



/ = is: and 
: as well 



where n — 1 for I — H , n — 2 
/ = i/2 n = 3 for / = HK and 
as / = G. These equations demonstrate the link be- 
tween the short distance expansion of the universal am- 
plitude function [Eq. (|4])] and the curvature expansion of 
the order parameter profile [Eq. pO|) ]. The temperature 
dependence of the short distance expansion coefficient 

functions ^±\r_i/S,±) follows from simple scaling con- 



siderations of the order parameter [31l . 13^ 
system under consideration: 



applied to the 



^i'^(r±/e±->0)-l + 4'i (^g) 



/ X 2/1/ / X D 

/ \ I/'' / \ 2/1/ 



(13) 



where I = p,H, K, H^, HK, H^, G. Since the last term 
in Eq. p3p represents the leading non-analytic contri- 
bution to the order parameter profile for t — > 0, one 
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has 



III. MEAN FIELD THEORY 

Within a field theoretical renormalization group ap- 
proach the scaling functions introduced above can be 
determined in lowest order perturbation theory within 
the framework of mean field theory corresponding to 
e = 4 — £' = 0. A renormalization group enhanced 
mean field theory is obtained by using for the scaling 
variables entering into scaling functions their full scaling 
form for e = 1. The standard Ginzburg-Landau fixed- 
point Hamiltonian for describing critical phenomena is 
given by [H,!!! 



(14) 



This Hamiltonian has to be supplemented by the bound- 
ary condition — +oo at the surface of the colloidal 
particle corresponding to the critical adsorption fixed 
point [35] . The integration runs over the volume acces- 
sible to the fluid and the parameter r is proportional 
to the reduced temperature t. Within mean field the- 
ory r = t/{i+f for T > and r = -\t\/{i^f/2 
for T < Tc with ^(f/CcT = V^. The coupling constant 
u > stabilizes the Hamiltonian 7i[0] for temperatures 
below the critical point (T < Tc) and (V(/))^ penalizes 
spatial variations of the order parameter. Within mean 
field theory fluctuations of the order parameter are ne- 
glected and only the configuration of the order param- 
eter with the largest statistical weight m = ^JuJ&{(j)) 
\/\T\/\t\/a is taken into account. After 



with y'u/6 

functional minimization one obtains the Euler-Lagrange 
equation 



Am — Tra + m 



(15) 



Equation psp can be solved numerically for arbitrary 
temperatures. For computational purposes it is conve- 
nient to choose a spheroidal coordinate system in which 
the surface of the ellipsoid corresponds to a constant 
value of one spheroidal coordinate (see, e.g., Ref. [36|). 
In the case of a generalized cylinder the Euler-Lagrange 
equation reduces to [13] 



9^ d-l d 

^-2-m(ri,i?,T) + , o TZ—'^i'^^'^^'^) 



dr 



rj_ + R dr_ 

R,T) + m^{r^,R,T), (16) 



where d — 1,2,3,4 denotes different types of gen- 
eralized cylinders. Using the short distance expan- 
sion [Eq. ^] for the universal amplitude function 
C+{r±/R) = r±m{r±,R,T = 0) [Eq. (0])] as input into 
Eq. (jl6p and equating terms of the same power in rj_/R 



leads to 

31 ^ 1 X 3 , , 

^^'=^216' ^^^'=2' ^^ = ^4' (^^^ 

where H = d-1, K^{d~ l)(d- 2)/2, and G = (d - 
l)(d-2)(d--3)/6. 



A. Curvature expansion coefficient functions 

The curvature expansion coefficient functions in 
Eq. (fTO|) are determined by solving Eq. (|16p for the four 
generalized cylinders with high symmetries and for arbi- 
trary values of r. In the following we restrict our presen- 
tation to the case T > T^. The first curvature expansion 
coefficient function is known analytically (see also 
the Appendix): 



(18) 



Our numerical data obtained from Eq. pop for d — 
2,3,4 can be used to determine individually all coefficient 
functions appearing in Eq. (jlOp. Moreover, the internal 
consistency of Eq. pH)) can be checked in the case of 



, (r^,i?),and0^^'^(r^,i?) 



the functions (/)). ' (r^ , R) 
in addition to the function (j}f\r±) for the half-space. 
This means that first each of these functions is deter- 
mined separately for each value of d, e.g., 0j d=2(''-L, R), 

ci=3 ('"-L I -R) ; 0j dLi{'''-L,R)- Thereafter these three 
functions as obtained this way are compared with each 
other. Internal consistency means that these three func- 
tions are identical within the numerical accuracy. For 
higher orders n > 3 in Eq. (|12p . there are contributions 
of more than two coefficient functions to the curvature ex- 
pansion. A consistency check for these functions would 
require the additional knowledge of the order parame- 
ter profiles around differently shaped surfaces. For com- 
parison we note that in the case of the aforementioned 
curvature expansion of density profiles of a hard sphere 
fiuid around an ellipsoid in 13 = 3 it was possible to 
check the internal consistency of the curvature expansion 
only in the case of the density profiles corresponding to 
4>\^\i'^,R) and (l)f\ri_) [13] ■ Hence the present calcu- 
lation in D = 4 provides a more stringent consistency 
test. Figures [5] (a) and (b) display the curvature expan- 
sion coefficient functions except for the well-known func- 
tion V+^(C+ = r^/£,+) for the half-space [Eq. One 
observes that all coefficient functions attain the contact 
value 1 at C+ = due to V+'(C+ ^ 0) - 1 = a{\Q+)^/'' 
in agreement with Eq. (|13p , where ^ = 1/2 within mean 

field theory and a^'* is a dimensionless coefficient (see the 
Appendix) . Moreover, all coefficient functions exhibit an 
exponential decay for large values of the scaling variable 
= rj^/^-|_. We have verified numerically the internal 
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FIG. 2: (a) [(b)] The curvature expansion coefficient functions 

ip+ \(+ = r±/^+) for / = H,K,H^ [HK,H^,G] of the order 
parameter profile for critical adsorption on curved surfaces 
according to Eqs. p2|l and (jlSp . The coefficient functions have 
been determined from the numerical solution of Eq. (|16p for 
d = 2, 3, 4 corresponding to three different types of generalized 
cylinders with high symmetries in spatial dimension D — 4. 
All coefficient functions attain the finite value 1 as ^ 
(see Eq. (US}). 



consistency of the first four terms of the curvature ex- 
pansion [Eq. pU])] as discussed above. 



B. Scaling functions 

We now turn our attention to the order parameter pro- 
files around the ellipsoidal particle shown in Fig. [TJ In 
Figs.[3](a) and (b) the scaling function = |r|~^/^m(T) 
of the order parameter profile is shown for two angles 
6 as obtained from Eq. (jisp together with the results 
of the curvature expansion [Eq. (fTO| ]. For 9 = 7r/2 
and R/^+ = 30 (lower curves in Fig. [3] (a)) the re- 
sults obtained from the curvature expansion are in agree- 
ment basically everywhere with the direct calculations 
of the scaling function. With increasing local curva- 
tures (see Fig. [ij, i.e., ioi 9 — and R/^+ — 30 in 
Fig. [3] (a) (upper curves) as well as for 6 = 0, 7r/2 and 




FIG. 3: (a) [(b)] The scaling function P+(C+ = 
r_L/£_+,R/i+,{e,},{Ss = Rs/R}) [Eq. i^] of the order pa- 
rameter profile for critical adsorption on an ellipsoid (see 
Fig. m with R/^+ = 30 [R/^+ = 15] and $2^63 = S4/2 = 1 
in dimension D — 4 for the angles 6 = 64, = ti/2 (lower 
curves) and 6 — 64—0 (upper curves). For symmetry rea- 
sons the order parameter profile is independent of the angles 
62 and 63. Direct numerical calculations (solid lines) accord- 
ing to Eq. (lisp are compared with the curvature expansion up 
to and including third order (dashed lines) as obtained from 
Eq. (|10p . The upper two curves have been shifted up by 1. 



— 15 in Fig. [3] (b), the third-order curvature ex- 
pansion (i.e., including terms up to n < 3) and the direct 
calculations deviate from each other for scaled distances 
C-H = ^ H,K,G because higher order terms in 

the curvature expansion are required. Nevertheless, the 
calculations support the use of the curvature expansion 
of the order parameter profile as a first approximation 
for large ellipsoids. 

The scaling function of the order parameter profile for 
an ellipsoid with two equal small {R = Ri = R2) and 
two equal large semi-axes (i?3 = R4) is shown in Fig. H) 
In the limit of a large aspect ratio i?4 ^ i? of the el- 
lipsoid the scaling function reduces to the one for a gen- 
eralized cylinder with d = 2 in Eq. (fT6|) . In the case 
that all semi-axes are the same the result for a sphere 
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FIG. 




r^/£,+,R/£,+,{es},{Ss = Rs/R}) (solid line) of the or- 
der parameter profile for critical adsorption on an ellipsoid 
with S2 — 5^/2 = 5^/2 = 1 and for the angles 6 = O4 = Tr/2 
(see Fig. [1} in dimension D — 4. For symmetry reasons 
the order parameter profile is independent of the angles 62 
and 63. The dashed and dotted lines represent the results 
for a sphere with equal semi-axes S2 = S3 = S4 = 1 and a 
generalized cylinder with 82 — 1 <^ 5^ — S4 corresponding 
to d = 4 and d = 2 in Eq. (|16p . respectively. The smallest 
semi-axis is fixed to R/^+ — 0.6 for the ellipsoid, the sphere, 
and the generalized cylinder. 




FIG. 5: The contour lines of the scaling function P+{(+ — 
r±/£,+,R/£,+, {63}, {5a = Rs/R}) of the order parameter pro- 
file for critical adsorption on an ellipsoid (thick line) in D = 4 
with R/^+ = 1 and S2 = S3 = 84 = 5. The center of the 
ellipsoid is located at the origin of the coordinate system and 
only the contour lines in the ("f — ("^ plane are shown, where 
= a;i/^+ and C^^ = 3^4/^+ (see Fig. [T}. The curves corre- 
spond to the following values of P+: 4, 1,0.1,0.01 (from the 
middle to the outside). 




10° 10' 10' 

A+1 

FIG. 6: The universal amplitude function C+(A = 
r±/R, {9s}, {Ss = Rs/R}) for an ellipsoid with 82 = S3 = 
S4/2 — 1 within mean-field approximation (i.e., D — 4) for 
the angles 6 = 64 = n /2 and 6 = Q (see Fig. [T|. For symme- 
try reasons the order parameter profile is independent of the 
angles 62 and 63. The short lines show the results obtained 
from the short distance expansion (up to and including the or- 
der indicated in Eq. Q). All curves share the common value 
C+(A = 0, {Os}, {Ss = Rs/R}) = c+ = y/2 corresponding to 
the universal amplitude for the half-space. 



is recovered. From Fig. [5] one can infer that the contour 
Hnes of the scahng function of the order parameter profile 
quickly exhibit spherical symmetry with increasing dis- 
tance from the surface of an ellipsoid. Along the radial 
direction 9 = (i.e., for = = 0) the scahng 

function decays faster away from the surface than along 
the radial direction = 7r/2 (i.e., for = Xi/£^+ = 0) 
because the local curvatures H , K, and G of the surface 
at = are larger than at 6' = 7r/2 (see Fig. [1]). For 
comparison we note that also in the case of generalized 
cylinders the decay rate of the scaling function increases 
upon increasing the local curvatures (see Figs. 5 and 6 in 
Ref. [131 )• 



C. Universal amplitude functions 

Whereas the preceding subsection has been focused 
on adsorption phenomena on ellipsoidal particles for 
T 7^ Tc, this subsection addresses the adsorption phe- 
nomena on ellipsoidal particles in the case that the sus- 
pending fluid is critical, i.e., at the critical composi- 
tion and at T = Tc- The universal amplitude func- 
tion C+(A = r±/R,{es},{Ss}) (see Eq. ([3])) is shown 
in Figs. El and [7] for two different ellipsoids. Accord- 
ing to Eq. (g]) the function C+(A, {6*4, {(5^}) starts at 
C+(A = 0,{9s},{Ss}) = c+, where c+ = within 
mean field theory. Figure [S] shows the numerical results 
as obtained from Eq. with C+ = r±m{T = 0) to- 
gether with the results according to the short distance 
expansion [Eq. (U)]. The calculations provide the range 
of validity A = r±/R < H,K,G of the short distance 



7 




lO" lO' lO' 10^ 

A+1 



FIG. 7: The universal amplitude function C+(A = 
r±/R, {6s}, {5a = Rs/R}) for two ellipsoids with two equal 
small and two equal large semi-axes ^2 = <53/25 = ^4/25 = 1 
(dotted lines) and S2 ~ Ss/A = 64/4 = 1 (dashed lines) in 
mean-field approximation (i.e., D = 4). The upper [lower] 
dotted and dashed lines represent the results for the angle 
9 — 71/2 [9 — 0] (see Fig. [!}. For symmetry reasons the 
order parameter profile is independent of the angles 62 and 
^3. The lower solid line represents the u niversa l amplitude 
function for a sphere of radius Reff = R\/S'l — 1/ Arcosh((54) 
with (54 = 4 as obtained from Eq. (|6]). The upper solid line 
follows from Eq. (|16p with d — 2 corresponding to a general- 
ized cylinder with radius R and two infinitely extended axes, 
i.e, 52 = 1, ^3 = ^4 = 00. All curves share the common value 
C+(A = 0,{9a},{Ss = Rs/R}) = V2. For large distances 
from the surface the universal amplitude function vanishes as 
A"'''^" with /3/i/ = 1 in the case of the sphere and the ellip- 
soids whereas for the generalized cylinder it tends to the finite 
value 1. 



expansion. For comparison we note that the curvature 
expansion for T ^ Tc provides a good approximation of 
the order parameter profile for all distances r± from the 
surface provided the local curvatures are smaller than the 
correlation length (see Fig. [3]). However, at the critical 
point the correlation length diverges and the curvature 
expansion, i.e., the short distance expansion, is valid only 
for distances from surface which are small compared with 
the local curvatures of the surface of the colloidal particle. 
Figure[7]demonstrates that the universal amplitude func- 
tion decays proportional to r ■^^^ with fijv = 1 for large 
distance rj^ from the surface. Far away from the surface 
of an ellipsoidal particle of arbitrary shape the univer- 
sal amplitude function is angle-independent and given by 
c'^^'^\rA_lReff) [Eq. ^] with an effective radius Reff 
that depends on the aspect ratios 5s of the ellipsoid (see 
the vanishing differences between the dashed lines and 
the lower solid line in Fig. [7| for large A). The expression 

Reff = R\/ S4 — 1/ Arcosh(54) for the effective radius fol- 
lows from a small particle operator expansion, analogous 
to the one presented in Refs. [U [13, HI] , to the system 
under consideration in this limiting case characterized by 
semi-axes of the ellipsoids which are much smaller than 
other lengths such as the correlation length or the dis- 



tance between the surface of the particle and the point 
at which the order parameter is monitored. 



IV. SUMMARY 

This work has been devoted to the investigation of crit- 
ical adsorption phenomena on ellipsoidal colloidal parti- 
cles [Fig. [Tj which are immersed in a fluid near critical- 
ity, i.e., t = {T - Te)/Te ^ in the case that the fluid 
is at the critical composition. The adsorption profiles 
are characterized by universal scaling functions P± for 
T ^ Tc [Eq. ([2])], involving the bulk correlations lengths 
and for T = hy universal amplitude functions C± 
[Eq. O]. 

For T — Tc the short distance expansion of the univer- 
sal amplitude function of the order parameter profile near 
a curved surface has been introduced [Eqs. ([3]) and (H])]. 
This expansion involves local curvatures of the surface 
[Eqs. ([5a|) ]. Also for temperatures T Tc a curvature 
expansion of the order parameter has been considered 
[Eq. pO)) ] which involves the local curvatures of the sur- 
face, too. The short distance behavior of the latter curva- 
ture expansion coefficient functions [Eqs. (fTT|) and (I12p ] 
is uniquely determined by a comparison of the curvature 
expansion [Eq. pO)l ] with the short distance expansion 
[Eq. (j4|)]. A numerical calculation within mean field the- 
ory [Eq. HH)] confirms this relation between the curva- 
ture expansion and the short distance expansion [Figs. [2] 
(a), (b) and [6]. 

For comparison we note that the contact values of cur- 
vature expansion coefficient functions of the density pro- 
file of a hard sphere fluid close to a hard convex wall 
are uniquely determined by exact statistical mechanical 
sum rules |39l . kol and the morphometric form of the 
grand potential [30, SH. If motion invariance, continuity, 
and additivity of the grand potential are satisfied, only 
four morphometric measures are needed to describe fully 
the influence of an arbitrarily shaped wall on thermo- 
dynamic properties of fluids consisting of spherical and 
non-spherical [l| particles. Hadwiger's theorem [H, IS] 
states that every motion-invariant, continuous, and ad- 
ditive functional in three dimensions can be expressed in 
terms of a linear combination of only four integrated geo- 
metric properties: the volume, the surface area, the mean 
integrated curvature, and the Euler characteristic of the 
wall. Thermodynamically away from bulk and surface 
phase transitions and for short-ranged interactions and 
correlations between the fluid particles, it is therefore ex- 
pected that the grand potential fulfills the requirements 
of Hadwiger's theorem. The value of the density proflle 
at a confining hard surface can be regarded as a ther- 
modynamic quantity because of exact sum rules imply- 
ing that the values of the curvature expansion coefficient 
functions of the density profile at a confining hard surface 
can be expressed in terms of thermodynamic properties 
such as the pressure of the bulk fluid, the fluid - planar 
wall surface tension, and the bending rigidity. These ar- 
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guments are expected to be not applicable to critical phe- 
nomena because intrinsic lengths in such systems reach 
macroscopic sizes so that the assumption of additivity is 
no longer valid. This has been confirmed by an explicit 
calculation of the excess adsorption for a sphere and a 
cylinder immersed in a fluid near criticality (see Figs. 8 
and 9 in Ref. [H). 

However, at T(. and close to a confining surface the 
curvature expansion coefficient functions of the order pa- 
rameter profile for critical adsorption are determined by 
the short distance expansion of the universal amplitude 
functions as discussed above; this expansion involves ar- 
bitrarily high orders. 

The curvature expansion for T ^ Tc provides a reli- 
able approximation of the order parameter profile for all 
distances from the surface provided the local curvatures 
are smaller than the inverse correlation length [Fig. [3], 
similar to the findings for the curvature expansion of the 
density profile of a hard sphere fluid close to a hard con- 
vex surface. However, at the critical point the correlation 
length diverges and the curvature expansion, which in 
this case reduces to the short distance expansion, is valid 
only for distances from surface which are small compared 
with the local radii of curvature of the surface of the col- 
loidal particle [Fig. [B]. In this sense at Tc the curvature 
expansion breaks down as a globally reliable approxima- 
tion. 

The scaling functions for critical adsorption on ellip- 
soidal particles of various shapes and sizes have been de- 
termined numerically for T ^ T,. within mean field theory 
[Figs. El m and[S]. In the case of very small or very large 
ratios of the semi-axes of the ellipsoid the results for a 
generalized cylinder such as a sphere or a rod are re- 
covered. At the surface the explicit form of the universal 
amplitude function C+ for an ellipsoid [Fig. [7] attains the 
value for the half-space and exhibits a power law limit- 
ing behavior far from the surface in accordance with the 
results for a sphere [Eq. ([6])] . 
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APPENDIX 

In this appendix we briefly discuss how the coefficients 
ipi'l- , '02^1 1 • • • tor I ~ p,H introduced in Eq. (fT^ can 
be calculated analytically within mean field theory, i.e., 
for D = A. Since the approaches for T > Tc and T < Tc 
are similar we restrict our presentation to the case T > 
Tc and drop the subscript '-)-' in order to simplify the 
notation. 

In the case of a generalized cylinder the scaling func- 
tion P{C = rj_/R,a = R/Cj = \T\-^l'^m{rA_,R,T) satis- 



fies the Euler-Lagrange equation [Eq. (|16p ] 

d-l 



C + a 



where the derivatives are taken with respect to the vari- 
able (. For a ^ 1 one may assume that the scaling 
function is analytic in a^^ such that it can be expanded 
into a Taylor series around a''^ = 0: 



F(C, a ^ (x.) = Po(C) + A(C)a~' + ' ' ' 



(A.2) 



By using the expansion {(+a)~'^ — a~^ + - ■ ■ in Eq. (jA.l[) 
in conjunction with Eq. (|A.2[) one obtains the familiar 
nonlinear differential equation for the profile near a pla- 
nar wall confining a semi-infinite system, 



with the solution 



^o(C) = - 



V2 



P 



(1) 



(A.3) 



(A.4) 



sinh(C) C 
and the linear differential equation 

Pi'iO + (d- l)P^(C) = Pi(C) + 3Po'(C)A(C) • (A.5) 

While the profile -Po(C) is well-known we now turn our 
attention to the calculation of Pi{C)- To this end we 
expand Pi{C) as 



PiiO^P, 



(0) 



-p^'\+p^^\'+ pf ^ e + Pt^ - 



(A.6) 



By inserting Eqs. (|A.4p and (|A.6p into Eq. (|A.5p and 
equating terms with the same power in C one derives the 
coefficients 



(0) 



3V2' ' 



(1) 



0, 



P 



(2) 



p(4) _ d-^ 



(A.7) 



6\/2 ' 72V2' 

(3) 

However, the value of P^ ' cannot be determined this 
way, i.e., for any value of P^^* Eq. (jA.SP is satisfied 
because of the known values of Pq^\ Plf\ Pq'^\ ■ ■ ■ . In 
Ref. [2^ the following analytic solution of Eq. (jA.Sp has 
been proposed: 



PiiO^id-1) 



cosh(C) — sinh(C) 



6V2sinh^(C) 
3Ccosh(2C) + 3sinh(2C) - 3Csinh(2C) 



2cosh(2C)-2-3C 
(A.8) 



This function can be expanded leading to the coefficients 
given m Eq. ([XT]) and to 



P 



(3) _ V2{d-1) 
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(A.9) 
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By inserting the scaling function P(C, a) from Eqs. (jA.2[) . 
(|A.6p . and (|A.7[) into Eq. ^ one obtains the expan- 
sion coefficient functions ipi^^ ? "02^'' 7 ' ■ ' for I = H in 
Eq. (fT3ll . We emphasize that a non- vanishing value of 
Pp"* would lead to an additional non-analytic term ~ t'^/^ 
in Eq. (fT3|) . However, according to the scaling consider- 
ations in Eq. the last term therein ~ |tp'' (= 



within mean field theory) represents the leading non- 
analytic contribution to the order parameter profile for 
t ^ 0. Therefore one concludes that = and the 
analytic expression in Eq. (jA.8[) does not capture fully 
the correct analytic properties of Pi {Q . Thus Eq. (jA.5|) 
has to be solved numerically. 
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